Aging of the Zero-Field-Cooled Magnetization in Ising Spin Glasses: 
Experiment and Numerical Simulation 
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A new protocol of the zero-field-cooled (ZFC) magnetization process is studied experimentally on an 
Ising spin-glass (SG) Feo.soMno.soTiOs and numerically on the Edwards- Anderson Ising SG model. 
Although the time scales differ very much between the experiment and the simulation, the behavior 
of the ZFC magnetization observed in the two systems can be interpreted by means of a common 
scaling expression based on the droplet picture. The results strongly suggest that the SG coherence 
length, or the mean size of droplet excitations, involved even in the experimental ZFC process, is 
about a hundred lattice distances or less. 
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PACS numbers: 75.10.Nr, 75.40. Gb, 75.50.Lk 

Spin glasses generally exhibit puzzling aging behav- 
ior such that their magnetic response in the spin-glass 
(SG) phase below its transition temperature T c depends 
crucially on their thermal history within the phase. For 
example, an isothermal zero-field-cooled (ZFC) magne- 
tization, mzFG,T(~ta, tw), is defined as the magnetization 
induced by a magnetic field h which is applied after a 
system is rapidly quenched from above T c to a tempera- 
ture T below T c and is aged for a waiting time of t w , and 
t a = t + t w with t being the time that the system ages 
at T under h. With increasing r, tozfc, T(ia 7 t w ) grows 
slowly. In contrast to ideal equilibrium situations, how- 
ever, it depends not only on r but also on the waiting 
time t w , namely, the larger t w , the slower the growth of 
TnzFC,T{t a ,t w ). This was the first observation of aging 
phenomena in spin glasses [Q. Since then most of the 
studies have focused on its t w -dependence. 

In the present work, on the other hand, we have stud- 
ied the ZFC magnetization focusing on its dependence 
on thermal processes after the field h is applied. We de- 
note it by mzFC,T[T](ta> tw) where t w symbolically repre- 
sents a common thermal history before h is applied, while 
T[t] does various heating processes after h is switched on, 
e.g., with a fixed heating rate but with an intermittent 
stop(s) at a different temperature (s). We have carried 
out this new protocol of the ZFC magnetization mea- 
surement both experimentally on an Ising SG material 
Fe .5oMno.5oTi0 3 §], hereafter referred to as FeMnTiO, 
and numerically on the 3D Gaussian Edwards- Anderson 
(EA) Ising SG model. Quite interestingly, we have found 
that rnzFC,T[T]{ta,tw) measured in the two systems are 
well interpreted by a unified way, or more explicitly, they 
are described by a common scaling expression which can 
be derived based on the droplet picture for the SG phase 
j| . The purpose of the present letter is to demonstrate 
this finding and to discuss its implications on the aging 
dynamics in spin glasses. 

According to the droplet picture, the mean size, R(t w ), 
of domains of an equilibrium state grows continuously 



when a spin glass relaxes with increasing t w . In fact 
its growth law in isothermal aging at a temperature T, 
denoted by Rt {t w Y, has been numerically studied by the 
various groups H-fl], and within the numerical accuracy 
all the results are consistent with each other. Here we 
quote our result M, 



R T (t w ) ~ (t w /r ) 



bfT/T c 



(1) 



where tq is the microscopic unit of time which is one 
Monte Carlo step (MCS) in the simulation, and b[ ~ 
0.17T C / J § ~ 0.153 with T c ~ 0.9J @], J being the vari- 
ance of the interactions. Furthermore we have shown j^] 
that the aging properties of the spin auto-correlation 
function are consistently explained when the same func- 
tional form is adopted for the maximum size of droplets 
which can be activated within a time scale of r, Lt(t), 



i.e., 



L t (t) ~ (t/to) 



b t T/T c 



(2) 



with the same constant bf . For an aging process in which 
the temperature is changed as in our new protocol, the 
maximum size of droplet excitations at time t a = r + t w , 
denoted by L(t), is evaluated by the successive use of 
eq.(||) as described by eqs.(||) and (||) below. 

In terms of L(t) introduced just above the common 
scaling expression we have found is written as 



l TnzFC.T[T}{ta,t w ) = xzFc(L(r),R(t w )), 



(3) 



where R(t w ) is the mean domain size at the instance 
when the field h is switched on. For isothermal pro- 
cesses under a field within the linear response regime the 
droplet theory || has already predicted a similar scaling 
form, i.e., tozfc, T(t a , t w )/h is a function of Rx{t w ) and 
Lt (t) . Equation (g) extends it to aging processes which 
involve temperature changes. Also, by the low temper- 
ature expansion of the droplet theory, mzFC,T{t a , t w )/h 
in the quasi-equilibrium regime with Rx(t w ) 3> Lt{t) is 
shown | |H not to depend explicitly on T. Equation (||) 
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shares the same property but for the whole range of r. 
In relation to this it is worth pointing out that the field- 
cooled (FC) magnetization, mFc(r), which is supposed 
to be limt a ^ 00 mzFC,T(^a, t w ), is nearly independent of T 
for the Ising spin glasses such as the FeMnTiO and the 
EA Ising model. 

The unified description of the ZFC magnetizations 
means that eqs.(|J) and (||) are quantitatively consistent 
with the measured results on the FeMnTiO spin glass 
when To ~ 10~ 12 s and the same value of bf as that of 
the EA-SG model are used. It is, to our knowledge, the 
first case in the study on aging in spin glasses that such 
a quantitative agreement is obtained in the direct com- 
parison of the experiment with the numerical simulation. 
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FIG. 1. The ZFC (o) and FC (x) magnetizations of the 
EA Ising SG model. The solid circles represent mzFc(T) with 
an intermittent stop of 10 4 MCS at T = 0.55. In the inset 
we show the corresponding results measured in the Ising SG 
Feo.soMno.soTiOs. The solid circles and squares represent 
tizfc(T) with an intermittent stop of 5 hours at 10. OK and 
15. OK, respectively. 

We have measured the induced magnetization of the 
single crystal FeMnTiO sample with T c =21. 5K @ by 
applying the field parallel to the hexagonal c-axis of the 
sample (see Q for further experimental details). The nu- 
merical simulations have been done on the 3D EA model 
with nearest neighbour interactions of mean zero and 
variance J which is the unit of energy. We have used 
the heat-bath dynamics and have simulated several hun- 
dred samples with a linear size L = 24 to make the error 
bars less than the size of the symbol in each plot. It is 
noted that the units of time and temperature used in the 
present work are respectively second and Kelvin in the 
experiment, and MCS and J, which we put unity below, 
in the simulation. Our processes measuring magnetiza- 
tions in the experiment [simulation] are as follows. The 
system was first cooled from 50K to 5K [2.5J to 0.3J] 
at a rate of 1.0 K/min [0.05 J/10MCS] under zero field. 
Once the lowest temperature was reached the field of 100 
Oe [0.2J] was applied and the sample was heated step- 
wise with a step of height 0.2K [0.01 J] and length 100s 



[100MCS]. The last 13s [10MCS] of the latter were used 
to measure mzFC,T[i-](*a,iw) which we regard as the ZFC 
magnetization at that temperature, tozfc(T). At 50K, 
in the experiment, we turned to decrease the tempera- 
ture and measured the FC magnetization, mpc(T), on 
cooling with the same rate and in the same field. In the 
simulation mFc(T) was measured on cooling the sample 
under h = 0.2J from 2.5J down to 0.3J with the same 
rate as for tozfc(T). 

The ZFC and FC magnetizations as well as the one 
with an intermittent stop of the real and model Ising spin 
glasses, shown in Fig. [j], look qualitatively quite similar 
to each other. Particularly, tofc(T) of both systems be- 
come nearly flat at T sufficiently lower than T c . We can 
also observe the following common feature in mzFc(T) 
with an intermittent stop. During the stop mzFc(T) in- 
creases, and when the heating is restarted, it continues to 
increase. However it increases significantly much slowly 
than mzFc(T) without the stop, and is caught up by the 
latter at a certain T below T c . 
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FIG. 2. Comparison of m Z Fc,T[i-](iaiiw) with one and two 
intermittent stops in the EA model (main frame) and the 
FeMnTiO (inset) [2]. Here they are denoted as mzFc(t) (see 
the text). 

In the inset [main frame] of Fig. |^ we demonstrate 
the most interesting observation in our new protocol. 
It is the comparison of the explicit time evolutions of 
mzFC, T[r]{ta,t w ) of the FeMnTiO spin glass [the EA SG 
model] with one intermittent stop at T = 15. OK [0.55J] 
and two stops at T = 14. 6K by 1 hour and at T — 15. OK 
[T = 0.50 J by 3000 MCS and at T = 0.55 J]. Here we 
rewrite them as mzFc(t) since we are interested in their 
time evolution. The origin of the time t is set at the 
instance when the system is heated up to 15. OK [0.55J] 
for the 1-stop process and to 14. 6K [0.50 J] for the 2-stop 
process. It is seen that, when the 1-stop curve (crosses) is 
shifted to the right by 2250 s [2100 MCS] (open circles), 
it almost perfectly collapses to the branch of T = 15. OK 
[0.55 J] of the 2-stop curve (solid circles). In other words, 
once rn Z Fc,T[T](ta,t w ) of two different thermal processes 
coincide with each other, their subsequent evolutions un- 
der a common condition of T and h are also identical. 
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Keeping in mind that these processes have a common 
R(t w ), or a common thermal history before h is applied, 
the results on both the real and model spin glasses are 
consistent with eq.(||) derived from the droplet picture. 
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FIG. 3. Evolution of the isothermal ZFC magnetization of 
the EA model at different temperatures but starting from a 
common Rr(tw) (— 2.1). In the inset, the same data are 
plotted against Lt{t) evaluated by eq.(2). 

To ascertain the scaling expression of eq.(||), and to 
demonstrate that the T-dependence of the scaling func- 
tion xzfc(L(t), R(t w )) comes out only through that of 
L(t), we have performed the simulation of isothermal 
ZFC magnetization processes at various temperatures 
but by choosing the waiting time f w after the instanta- 
neous quench such that Rx(t w ) given by eq.(|l|) become 
the same for all temperatures. The results are shown in 
Fig. |^ where mzFC ,t (t a ,t w ) are plotted against the time r 
after the field h = 0.2 is switched on. The chosen Rx(t w ) 
is ~ 2.1 and the corresponding t w at each temperature 
are indicated in the figure. 

The evolutions of mzFC.T{t a , t w ) in Fig- ^ depend on T 
except at r ~ 1. When, however, the same data are plot- 
ted against Lt{t) evaluated by eq.(||), they are nicely 
superposed upon each other as shown in the inset of 
Fig. ^. This result strongly supports the scaling expres- 
sion of eq.(||) in the following two respects: the thermal 
history before the field is applied is in fact renormalized 
to R(t w ), and the size Lt(t) of droplet excitations which 
are dominant in the time scale r after the field is ap- 
plied is properly given by eq.(||). We also note that, 
within the accuracy of the present numerical analysis, 
the scaling behavior is ascertained not only in the quasi- 
cquilibrium regime, where Lt{t) < Rr(t w ) or t <C t w , as 
predicted by the droplet theory ||U , but also in the out- 
of-equilibrium, or aging regime, where Lj>(t) > i?r(^w) 
or t > f w . In the latter regime Lt{t) plays a role of the 
mean domain size at time t a in place of RT{t w ). 

Lastly let us extend the scaling analysis done above to 
our new experimental protocol in which the temperature 
is changed. For this purpose it is important to note that, 
so long as the system is in the SG phase, the mean do- 



main size R(t w ) of the EA model analyzed numerically 
ever grows continuously even when the temperature is 
changed discontinuously |ll| . It is then natural for us to 
expect that this is also the case for L{t) under a fixed 
field h, so that we may successively use eq.@ as follows. 
Suppose that temperature is changed as T t — Tj_i + AT* 
and a period of time d,Ti is spent at each temperature Ti. 
We assume that the length scale achieved after the step 
at Ti is given as 



ft 



off 



dn 



\b t Ti/T c 



(4) 



where t° is an effective elapsed time to take into account 
the thermal history up to the instance just before the 
temperature is changed to Ti. It is the time needed to 
achieve the length scale L^-i reached after the [i — l)-th 
step by isothermal aging at 



cffN&fTi/i; 



(rf) 



Li-i 



(5) 



Combining these we obtain a recursion formula for Li 
and/or r? ff which is solved for any given temperature 
process with L(t — 0) = 0. 
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. . j. 4. Scaling plot of jtizfc t[t] (tavfw) of the EA model 
against L(r) evaluated by eqs.(Q) and (|^). The solid line rep- 
resents the one without intermittent stop, while the symbols 
are those at the stop at various temperatures T where r cff 
is the effective elapsed time at the instance when the system 
just reaches T. The data at T = 0.55 with r cft = 2100 MCS 
represent the second stop in the 2-stop process. The dashed 
line is one of the curves shown in the inset of Fig. 3. 

We have first applied the above-mentioned analysis to 
TO ZFC,T[T](ia ; t w ) of the EA model in the processes with 
T-changes shown in Figs. [I] and ^. The results are shown 
in Fig. y. All mzFC,T[T] {ta, tw) nicely lie on a univer- 
sal curve, indicating that the scaling of eq.(^), combined 
with the above analysis on L(r), does work well. In the 
figure we also draw one of the curves in the isothermal 
process shown in the inset of Fig. [| It has the differ- 
ent thermal history before the field is applied, i.e., the 
different R(t w ), from those of the T-change process. It 
significantly deviates from the universal curve for the lat- 
ter as expected from eq.(H). 



3 



Surprisingly, the rnzFC,T[T](t a ,t w ) curves measured in 
the FeMnTiO spin glass and shown in the inset of Figs, [j] 
and H obey the same scaling law as those of the EA model. 
Indeed, their scaling plots demonstrated in Fig. || are ob- 
tained by making use of eq.(^) with the identical expo- 
nent of bfT/T c but with r = 10 -9 s and 10 -12 s, instead 
of tq = 1 MCS for the EA model. The two portions with 
the dense data points at around rn ZFC T ^(t a , t w ) ~ 0.14 
and 0.20 represent the results with intermittent stops 
at 9.6K, 10.0K and at 14.6K, 15. OK, respectively. At 
a glance we see that the scaling with To = 10~ 12 s is 
the better. We have not performed further analysis on 
the sensitivity of the collapse of the curves to the pa- 
rameter values, because the evolution law of eq.(||) or 
(Q) already involves certain uncertainty. For example, 
Rr(t w ) extracted by the ac susceptibility in a sufficiently 
wide frequency range is better expressed by a function of 
\nt w ]l2] |. Furthermore we are faced to similar but more 
serious problems in the temperature range close to T c , 
such as the crossover to the critical regime ||Jl^] , whose 
analysis has been omitted in the present work. How- 
ever we consider that the scaling results obtained above 
are enough for us to claim that the droplet picture ex- 
pressed by eq. (|^) is consistent with the aging behavior of 
the ZFC magnetization in the FeMnTiO Ising spin glass, 
and that the length scale involved in the process is about 
a hundred lattice distances or less. They also imply the 
continuous growth of L(t) and so R(t w ), without suffer- 
ing from the chaotic effect at least in the present 
experimental protocol. 
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FIG. 5. Scaling plots of fn ZF c,T[T](*a,£w) of the FeMnTiO 
spin glass shown in the inset of Figs. 1 and 2. The abscissa 
L(t) is evaluated by eq.(2) with to = 10 -9 s and 10 -12 s. 

Another interesting result of the present work is that 
the evolution law of Lt{t) given by eq.(^), or the growth 
law of i?r(tw) given by eq.([j]), almost coincides with the 
growth law of the SG coherence length, £,r(t), which 
Joh et al. |l^,0] has extracted from the thermoremanent 
magnetization(TRM) of the spin glasses CuMn, AgMn 
and Cdlnp^Cr! 7S4, though none of them is an Ising spin 
glass. Their results imply that the evolution law of Lt{t) 
depends on the field at least in the the crossover region 



between the quasi-equilibrium and the aging regimes. In 
the present work, however, such a field effect has been 
neglected regarding it as a higher order effect. 

To conclude we have demonstrated that the ZFC mag- 
netizations measured in both real and model Ising spin 
glasses are described by a common scaling expression 
which is consistent with the droplet picture on the SG 
aging dynamics. The result indicates that the numeri- 
cal simulation on the model spin glass supplies us not 
only qualitative but also quantitative information on ag- 
ing phenomena in real spin glasses. 
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